A geospatial analysis of energy transition mining in the U.S.

mines <- minesRaw %>% #filter to key energy transition minerals
  filter (CRITICAL_MINERAL %in% c("Aluminum",
                                  "Cadmium",
                                  "Cobalt",
                                  "Copper",
                                  "Dysprosium",
                                  "Gallium",
                                  "Indium",
                                  "Lithium",
                                  "Manganese",
                                  "Neodymium",
                                  "Nickel",
                                  "Silver",
                                  "Selenium",
                                  "Tellurium",
                                  "Rare-Earth Elements"))

minesKeyMinerals <- minesRaw %>% #filter to key energy transition minerals
  filter (CRITICAL_MINERAL %in% c("Cobalt",
                                  "Lithium",
                                  "Tellurium",
                                  "Rare-Earth Elements"))

map3 <- ggplot() + 
  geom_polygon(data=state, #add base map
               color = "black", 
               fill="white",
               aes(x=long, y=lat, group=group)) + 
  geom_point(data=minesKeyMinerals, #add mines
             aes(x=LONGITUDE, y=LATITUDE, color=CRITICAL_MINERAL,
                 text=paste("",CRITICAL_MINERAL))) +
  theme_minimal()+ #remove grey background
  theme(legend.position="none", #remove legend, axes labels/ticks
        axis.title.x=element_blank(), 
        axis.text.x=element_blank(), 
        axis.ticks.x=element_blank(),
        axis.title.y=element_blank(), 
        axis.text.y=element_blank(), 
        axis.ticks.y=element_blank()) + 
  ggtitle('"Critical Mineral" Mines & Deposits in U.S.:<br>Cobalt, Lithium, Tellurium, Rare Earths') +
  coord_fixed(1.3) #adjust aspect ratio

ggplotly(map3, tooltip="text") %>% #make interactive
  style(hoverinfo="none", traces=1)

# map4 <- ggplot() + 
#   geom_sf(data=county_simpl, #add base map
#                color = NA, 
#                aes(fill=estimate)) + 
#   coord_sf()+
#   scale_fill_viridis_c(option = "magma") +
#   geom_point(data=minesKeyMinerals, #add mines
#              aes(x=LONGITUDE, y=LATITUDE, color=CRITICAL_MINERAL,
#                  text=paste("",CRITICAL_MINERAL))) +
#   theme_minimal()+ #remove grey background
#   #theme(legend.position="none", #remove legend, axes labels/ticks
#         # axis.title.x=element_blank(), 
#         # axis.text.x=element_blank(), 
#         # axis.ticks.x=element_blank(),
#         # axis.title.y=element_blank(), 
#         # axis.text.y=element_blank(), 
#         # axis.ticks.y=element_blank()) + 
#   ggtitle('"Critical Mineral" Mines & Deposits in U.S.:<br>Cobalt, Lithium, Tellurium, Rare Earths') 
#  # coord_fixed(1.3) #adjust aspect ratio
# 
# ggplotly(map4, tooltip="text") %>% #make interactive
#   style(hoverinfo="none", traces=1)

#buf35miles <- st_buffer(iowa_methodist, dist = 50)